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Abstract 

In this note, starting with a httle-known result of Kuo, I derive a recurrence 
relation for the Bernoulli numbers B2n, n being any positive integer. This new 
recurrence seems advantageous in comparison to other known formulae since 
it allows the computation of both B^n and -B4n+2 from only Bq, B2, . . . , -B2n- 
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The Bernoulli numbers n being a nonnegative integer, can be defined 
by the generating function 



e^ — 1 n\ 



ri=0 



The first few values are well-known: Bq = 1, Bi = — |, and B2 = \- It is 
also well-known that = for odd values of n, n > 1. For even values 
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of n (= 2m), n > 1, the numbers form a subsequence of non-null real 
numbers such that (—1)™+^ > 0. In other words, the entire sequence of 
Bernoulli numbers up to any i?„, n > 1, consists of Bq and given above, 
and the preceding numbers B2mi m = 1, . . . , [n/2\. The basic properties of 
these numbers can be found in Sec. 9.61 of 

The numbers Bn appear in many instances in pure and applied mathe- 
matics, most notably in number theory, finite differences calculations, and 
asymptotic analysis. Therefore, the efficient computation of the numbers Bn 
is of great interest. To this end, recurrence formulae were soon recognized 
as the most efficient tool One of the simplest such relations is found by 
multiplying both sides of Eq. ([1]) by e^ — 1, using the Cauchy product with 
the Maclaurin series for e"^ — 1, and equating the coefficients of the powers 
of X, which results in 

j=o ^ ^ 

This kind of recurrence formula has the disadvantage of demanding the pre- 
vious knowledge of all Bq, Bi, . . . , Bn-i for the computation of Bn- In search- 
ing for more efficient formulae, shortened recurrence relations of two different 
types have been discovered. The first type consists of the so-called lacunary 
recurrence relations, in which Bn is determined only from every second, or ev- 
ery third, etc., preceding Bernoulli numbers (see, e.g., the lacunary formula 
by Ramanujan [3;]). The second type demands the knowledge of only the 
second-half of the Bernoulli numbers up to Bn~i in order to compute Bn 4|. 



For an extensive study o 
Bernoulli numbers, see (5|. 



these and other recurrence relations involving the 
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Here in this note, I apply the Euler's formula relating the even zeta value 
C{2n) to i?2n to a little-known recurrence formula for C{2n) obtained by 
Kuo [6], in order to convert it into a recurrence formula for i?2n- By doing 
this, in fact I introduce a third type of recurrence relation for the Bernoulli 
numbers, in the sense that it allows us to compute both i?4„ and -B4„+2 from 
only the first-half of the preceding numbers, i.e. Bq, B2, . . . , B2n- Then, the 
efficiency of this new recurrence formula certainly surpasses that of most 
known formulae, mainly for large values of n. 

For real values of s, s > 1, the Riemann zeta function is defined as ({s) : = 
Yl'^=i ^/n'"^- In this domain, the convergence of this series is guaranteed by 
the integral test l] For positive even values of s, one has the well-known 
Euler's formula (1740) [Tj: 

o2n— 1 f) 

Ci2n) = i-ir-' (3) 

which, in face of the rationality of every Bernoulli number Bn, yields C{2n) 
as a rational multiple of vr^". The function C{s), as defined above, can be 
extended to the entire complex plane (except the only simple pole at s = 1) 
by analytic continuation, which yields ((0) = — The reader should note 
that Eq. ([3]) remains valid for n = 0. 

These are the necessary ingredients to state our first lemma, which comes 



^For s = 1, one has the harmonic series X^^^i which diverges to infinity. 
^Alternatively, we can use the globally convergent series 

1371^ E,T=0 2^ ELo(-i)Mfc) (fcTTF' *° ^''^'^ (1^^^) H' ^^'""^ ^^^^^ 

for all complex numbers s ^ 1 + i: being any integer. At s = 0, it reduces to 
C(0) = -EZol^ELoi-^ni) - -^-^Er=i2^ELo(-l)M^)- Of course, 
ELo (-1)' (fe) = for all n > 0, so C(0) - -i. 
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from a little-known recurrence formula by Kuo (1949) 6|, just written directly 
in terms of the Riemann zeta function. 



Lemma 1 (Kuo's recurrence formula for ({2n)). For any positive in- 
teger n, one has 



i2n-l7r2" 



\n/2\ 



C(2n) 



4(n- l)!'(2r2- 1) (n - 1)! ^ ^ ' 



k=0 



{n - 2k)\ {2n - 2k) 



[n/2\ ln/2\ 



(27r 



fc=0 j=0 



> - 2k)\ {n - 2j)! (2n - 2A; - 2j + 1) 



This recurrence formula is proved in [6| by developing successive inte- 
grations (from to x) of the Fourier series sin (?2x)/n = (vr — x)/2, 
which converges for all positive real x < 2 vr, and then applying the Parseval's 
theorem. 

We are now in a position to prove the following theorem. 



Theorem 1 (Recurrence relation for B2n)- For any positive integer n, 
one has 



B2n = {-ir-' a 

[n/2] „ [n/2\ 



ln/2] 



B2k 



k=0 



{2k)l {n - 2ky. {n - k) 



+ (2nV V — V ^ - 

^ {2k)\{n-2k)\ (2j)!(n-2j)! 2n - 2A; - 2j + 1 



(4) 



where ci„ — ^ |„_]^)i2 and — 2(n—i)\ ' 



(n-iy. 
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Proof. By dividing both sides of the Kuo's recurrence formula, as given in 
Lemma [H by tt^"", one has 



C(2n) 



i2n-l7r^ 



TT 



2n 



+ 



1 



A{n-l)f{2n-l) {n - 1)1 



ln/2] 

E (-1)' 

k=0 



-)2n-2k 



ln/2\ \n/2\ 



C(2fc) 

7r2fc {n - 2k)\ {2n - 2k) 

22n-2k-2j+l 



k=0 j=0 

From Euler's equation, Eq. ([3]), one knows that 



7r2i (n - 2k)\ {n - 2j)! (2n - 2A; - 2j + 1) 



• (5) 



-ly 



-1 2^"'-! B2„ 



(2m)! ' 

which is vahd for all integer m, m > 0. By substituting this in Eq. ([5]), one 
finds, after some algebra. 



B. 



2n 



-1) 



ri-l 



2A; 



1 



i2n)\ A{n- l)\' {2n - 1) {n - 1)! ^ (2fc)! {n - 2k)\ {2n - 2k) 

\n/2\ \n/2\ 



+(-1)"-' E E 



B2k 



1 



(2A;)! (2j)! (n - 2^)! (n - 2j)\ {2n -2k- 2j + 1) 



■ (6) 



k=0 j=0 

Now, put (— 1)"""^ in evidence and multiply both sides by {2n)\. This yields 



B 



2n 



-1) 



n-1 



{2n) 



(2ra)! ^ B2k 



A{n- l)f {2n - 1) (n - 1)! ^ (2A;)! - 2k)\ {2n - 2k) 



\n/2\ \n/2\ 

+i2n)\ y: e 

fc=0 j=0 

which readily simplifies to 



B2k B2j 



{2k)\ {2j)\ {n - 2k)\ {n - 2j)\ {2n -2k- 2j + 1) 



, (7) 



B 



2n 



n-1 



-1) 



\n/2\ 



n{2n-2)\ {2n)\ 



[n/2\ 

E 



2k 



2{n- 2{n- 1)! ^ (2A;)! (n - 2k)\ {n - k) 



+ {2n)\ J2 



B 



2k 



[n/2i 

E 



B. 



2i 



(2A;)! (n - 2k)\ ^ {2j)\ {n - 2j)\ 2n - 2k - 2j + 1 



• (8) 



□ 
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In fact, the auxiliary terms a„ and 6„ in our Theorem [T] can be cast in a 
more suitable form for computational purposes. For a„, one has ai = | and 

fi (2n — 2)' n -A- 

for all n > 1. For 6„, one has 61 = 1 and 

(2n)(2n-l)- ... -n (n-1)! 
bn = ^ '- = n'{2n-l). ....{n + l) 

Y[{2n-m), (10) 



n-l 

,2 



= n 
for all n > 1. 

Since Bq = 1 and the recurrence relation in Theorem [1] has only the four 
basic numeric operations (+, — , x, and -=-), it is straightforward to show, by 
induction on n, that every B2n, n > 1, is a rational number, though this is a 
well-known characteristic of these numbers (see, e.g., Ref. jsl). 
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